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ABSTRACT 

Observations of the microwave sky using the Python telescope in its fifth season of operation at the 
Qv| . Amundsen-Scott South Pole Station in Antarctica are presented. The system consists of a 0.75 m off- 

axis telescope instrumented with a HEMT amplifier-based radiometer having continuum sensitivity from 
37-45 GHz in two frequency bands. With a 0.91° x 1.02° beam the instrument fully sampled 598 deg^ 
iL^i of sky, including fields measured during the previous four seasons of Python observations. Interpreting 

the observed fiuctuations as anisotropy in the cosmic microwave background, we place constraints on 
the angular power spectrum of fluctuations in eight multipole bands up to Z ~ 260. The observed 
spectrum is consistent with both the COBE experiment and previous Python results. There is no 
significant contamination from known foregrounds. The results show a discernible rise in the angular 
^S| ■ power spectrum from large {I ~ 40) to small [I ^ 200) angular scales. The shape of the observed power 

^ ' spectrum is not a simple linear rise but has a sharply increasing slope starting at I ~ 150. 

Subject headings: cosmic microwave background - cosmology: observations 
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two feeds observe the same elevation and are separated by 
2.80° on the sky. These beams are scanned horizontahy 
across the sky by a large rotating vertical flat mirror, the 
chopper, at 5.1 Hz. The new scan strategy motivated two 
changes from the instrument configuration described in 
Kovac et al. (1999): the frequency response of the data 
system was extended by switching to 100 Hz-roUoff an- 
tialiasing Bessel filters, and the data recording rate was 
correspondingly increased, to 652.8 samples/sec for each 
channel. 

3. CALIBRATION 

As in previous Python seasons, the primary DC calibra- 
tion of the detectors was derived using liquid nitrogen, liq- 
uid oxygen, and ambient temperature thermal loads exter- 
nal to the receiver (Dragovan et al. 1994, Ruhl et al. 1995). 
Load calibrations were performed approximately once per 
day, and gains were found to be consistent over the entire 
season to within ± 2%, with no discernible trends. Gain 
compression, which was a source of systematic uncertainty 
in the calibration of Pyl-Pylll, is measured to be negligible 
for the Python HEMT receiver. Systematic uncertainty in 
the DC load calibration is estimated to be ± 10%. 

Several efficiencies must be estimated to relate the load 
calibrations to celestial response in the main beam, which 
account for power losses in the atmosphere, in the side- 
lobes, and in the telescope, and they are calculated using 
data from skydips and from various beam measurements. 
The resulting systematic calibration uncertainty of ^w 
is asymmetric, due to the fact that the individual losses 
are small positive numbers and hence the errors in their 
estimation follow skewed distributions (Kovac et al. 1999). 

The dynamic response of the system was calculated from 
laboratory measurements of the transfer functions for the 
AC coupling and antialiasing filters in the data system, 
and confirmed on the telescope by comparison of obser- 
vations made of the moon using normal and slow chop- 
per speeds. The response speed of the detectors is not 
a concern for this calibration or for its uncertainty An 
appropriate response correction factor is applied to each 
modulation of the data. The uncertainty on these factors 
is small, and is dominated by a ± 5% systematic uncer- 
tainty on their common normalization. 

The overall uncertainty in the calibration of this dataset 
is estimated to be ^^jr- Antenna temperature has been 
converted to units of STqmb throughout. 

4. OBSERVATIONS 

Two regions of sky were observed: the PyV main field, 
a 7.5° X 67.7° region of sky centered at a = 23.18'', 6 — 
—48.58° (J2000) which includes fields measured during the 
previous four seasons of Python observations and a 3.0° x 
30.0° region of sky centered at a = 3.00'', S = -62.01° 
(J2000), which encompasses the region observed with the 
ACME telescope (Gundersen et al. 1995). The total sky 
coverage for the PyV regions is 598 deg^, greater sky cov- 
erage than previous degree-scale CMB experiments. The 
combined absolute and relative pointing uncertainty is es- 
timated to be 0.15°, as determined by measurements of the 
moon and the Carinae nebula {a — 10.73h,S ~ —59.65°). 
The PyV beam is well approximated by an asymmetric 

^°'°^ degrees (az x el). 



The beam is determined from scans of the Carinae nebula 
and the Moon. Both PyV regions are fully sampled with 
a grid spacing of 0.92° in elevation and 2.5° in right as- 
cension, corresponding to a distance of 1.6° on the sky at 
a declination of —50°. The telescope is positioned on one 
of the fields and the chopper smoothly scans the beams in 
azimuth in a nearly triangular wave pattern. The chop- 
per throw is 17° in azimuth, corresponding to 11° on the 
sky at a declination of —50°. A total of 309 fields are 
observed in 31 sets of 5-17 fields. Some of the fields are 
observed in more than one set. There are 128 data sam- 
ples for each detector channel in a complete chopper cycle, 
and 164 chopper cycles of data are taken of a given field 
before the telescope is positioned on the next field in the 
set. One data file consists of 164 chopper cycles for each 
field in the set. Approximately 10 hours of good data (100 
files) are taken of a set of fields before the telescope moves 
on to the next set of fields. 

5. DATA REDUCTION 

After cutting 45% of the data for weather and tracking 
errors, 389 hours of data remained for use in the CMB 
analysis. The data are modulated using 



Mm{9) = cos{Tmre/ec) x 
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H{0) m 
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where m is the modulation number, 9 is the chop- 
per angle, 9c is the extent of the chopper throw and 
H{e) = 0.5(1 - cos{2ti9/9c)) is a Hann window. (Fig. 1). 
The TO = 2 ... 8 modulations are apodized with the Hann 
window in order to reduce the ringing of the window func- 
tions in ^-space. Data taken during the right- and left- 
going portions of the chopper cycle are modulated sepa- 
rately, to allow for cross-checks of the data. Sine modu- 
lations are not used in the analysis because they are anti- 
symmetric and are thus sensitive to gradients on the sky. 
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Fig. 1. — The data are modulated with cosines which have been 
apodized by a Hann window. Three of the modulations are shown. 

Data in a given file, field, channel, and modulation are 
co-added over all chopper cycles. A chopper synchronous 
offset, due to differential spillover past the chopper, is re- 
moved from each data file by subtracting the average of 
all of the fields in a file. This is not just a DC offset; there 
is an offset removed for each modulation (or equivalently, 
for each sample of the chopper waveform) and channel. 
This subtraction must be accounted for by adding a term 
to the covariance matrix (see section 6). When the mod- 
ulated data are binned in azimuth, a periodic signal due 



to the 12 panels of the ground shield is evident, especially 
on larger angular scales. The signal of period 30° is fit for 
an amplitude and is subtracted. Removal of the ground 
shield offset has less than 4% effect on the final angular 
power spectrum (section 6) in all modulations. 

After the data have been modulated and offsets re- 
moved, the right and left-going data, which have been 
properly phased, are co-added, as are channels which ob- 
serve the same points on the sky. Since the two feeds ob- 
serve different points on the sky they cannot be co-added; 
the theoretical and noise covariances between them are 
included in the likelihood analysis of the angular power 
spectrum. 

The South Pole Station power nominally operates at a 
frequency of 60 Hz, so there is possible contamination at 
60 Hz and its harmonics. However, the chopper runs at 
5.1 Hz, which is incommensurate with 60 Hz. 

Our noise model assumes the covariance between fields 
taken with different sets of files is negligible because of the 
chopper offset removal and because of the long time (at 
least 10 hours) between measurements. An analysis com- 
paring the noise covariance estimated from data which had 
not yet been co-added over all cycles and the noise covari- 
ance estimated from data which had been co-added over 
all cycles indicates that PyV noise is dominated by de- 
tector noise. However, the long term drifts due to the 
atmosphere, which add to the variance as well as induce 
small correlations between fields taken with the same set 
of files, are important, especially for power spectrum es- 
timation (section 6). The noise covariance between fields 
taken with the same set of files is first estimated by taking 
the usual covariance on the co-added data, but because 
there were typically only 100 files taken for each field, the 
sample variance on the noise estimate is ~ (100)^^^, or 
10%, which will severely bias estimates of band power. To 
obtain a better estimate of the noise, we averaged the vari- 
ances for each set of files and then scaled the off-diagonal 
elements of the covariance to the average variance in a 
given set based on a model derived from the entire PyV 
data set. 

Several self-consistency checks were performed on the 
data set, using our best estimate of the noise covari- 
ance. The x^ = d*C~^d, where d is the data vector and 
C = C'^ +C^ +C'^ is the total covariance (see section 6), is 
consistent with its expected value (the number of degrees 
of freedom) in all modulations. The value of the proba- 
bility enhancement factor, /3, (Knox et al. 1998) between 
data from the two feeds falls within the expected range 
for all modulations. Finally, the data set was transformed 
into the signal-to-noise eigenmode basis. In that basis, 
C^^/^d should be Gaussian distributed with a = 1 and 
a total area equal to the number of degrees of freedom. 
Histograms of C~^/^d are consistent with such a Gaus- 
sian distribution for all of the modulations. We performed 
these same tests using the preliminary noise covariance, 
which suffers from sample variance, and found that the 
consistency tests failed. The tests indicate that the data 
set used is self consistent and that our best estimate of the 
noise covariance is indeed a good model for the noise. 

6. DATA ANALYSIS AND RESULTS 

CMB angular power is usually expressed in terms of 
angular multipoles, C;. A flat power spectrum is one for 



which C = (^(Z-f l)C;/27r) is constant. For each of the eight 
modulations, we compute the likelihood (£) as a function 
of C. The theoretical covariance matrix, C'^ , needed to 
compute C depends on the amplitude of C and the exper- 
imental window functions. The window functions are a 
measure of experimental sensitivity as a function of angu- 
lar scale I (Fig. 2). They are generated from the exper- 
imental beam map, modulations, and observing strategy 
and are given by 



Wu, = (1/2^) 



f d(f>e-'^<'''-''^'>B{k)B*{k) (2) 



where -B(k) is the Fourier transform of the beam map 
for the given modulation, Xj the position of field i, and 
k — Z(cos(/>, sin^). These functions are computed for all 
pairs of fields and channels. 




Fig. 2. — Diagonal window functions. The unapodized cosine 
modulation is plotted with a dashed line and the apodized cosine 
modulations are plotted with solid lines. 

The subtraction of the chopper synchronous offset can 
be accounted for with an additional term in the covariance 



matrix (e.g., (Bond et al. 1998)) and in this case is given 
by C*^ = {K/Nfiies) for fields taken with the same set of 
files and zero for fields not taken with the same set of files. 
iVfiios is the number of files which observed the fields. Tak- 
ing K to be large ensures that we have no sensitivity to 
modes of the data that could have come from the chopper. 
The CMB power spectrum is shown in Fig. 3 and is 
given in Table 1. The band powers are calculated for 
each modulation separately. This neglects the correla- 
tions between modulations. Work is in progress to fit 
the amplitudes in a series of bands to all of the modu- 
lations simultaneously. The calibration uncertainty allows 
all band powers to shift by the same amount (i.e. the cal- 
ibration errors are correlated). Given the uncertainty in 
the beam size of approximately 0.015°, the band power 
for a given modulation can move roughly by a factor of 
exp(±Z(0.425)(0.015)(7r/180)), only a 3% effect at / = 200. 
Flat band powers calculated from the subset of PyV data 
in the Pylll region of sky are consistent with the PylH 
flat band powers to within the uncertainties. Given the 
^ifjy calibration uncertainty and the 3/iK statistical un- 
certainty, the lowest I PyV modulation agrees with the the 
smallest scale COBE measurements. 



Table 1 

Results of likelihood analysis. Band powers are in iJ,K. 
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7. FOREGROUNDS 

The PyV data are cross-correlated with several fore- 
ground templates in order to set limits on possible fore- 
ground contamination. The templates used are the 
Schlegel et al. (1998) 100 micron dust map, which is based 
on IRAS and DIRBE maps, the Haslam et al. (1974) 408 
MHz survey (synchrotron), and the PMN survey (point 
sources). Each foreground template map is smoothed to 
PyV resolution, pixelized and modulated according to the 
PyV observation scheme. 

Two templates are created for the PMN survey. We call 
one PMN, which is converted to 5Tcmb using the spec- 
tral indices given in the survey. The other we call PMNO, 
which is converted to a flux at 40 GHz assuming a flat 
spectrum extrapolated from the flux measurement at 4.85 
GHz. The assumption of a flat spectrum is conservative in 
that it is likely to over-estimate the flux at 40 GHz. Nei- 
ther case is correct, since spectral indices have not been 
measured for all of the sources, in which case a flat spec- 
trum is assumed, but we do know that some of them are 
not flat, so a flat spectrum will be inappropriate. The two 
cases cover a reasonable range of possibilities. 

For each modulation and foreground template, a corre- 
lation coefhcient and uncertainty are calculated following 
de Oliveira-Costa et al. (1997). A weighted mean and un- 
certainty over all of the modulations for each foreground 
are given in Table 2. In all cases there is no clear detec- 
tion of foregrounds. The RMS of each modulation of each 
foreground was calculated and then multiplied by the cor- 
responding 1(7 error bar in Table 2 in order to estimate an 
upper limit on the foreground contribution to CMB band 
power. The limits on contributions from foregrounds are 
given in Table 3 and are at least ^ lOx smaller than the 
measured CMB bandpowers. 



Table 2 

Correlation coefficients and uncertainties for weighted means. 



Dust 
^iKiMJy/sr)' 
^3"^ 



18 



Haslam 

rXe 0.012 ± 0.024 



-2.0 



PMN 



PMNO 
fiK{MJy/sr)~ 



195 ± 385 



Table 3 

Upper limits on foreground contribution. All units are fiK . 
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If the diffuse morphology of the sky is not constant as 
a function of wavelength, then these templates do not re- 
veal all of the foreground contamination and more could 
be hidden in the PyV data. A combined analysis of Pylll 
at 90 GHz and PyV at 40 GHz would constrain the fore- 
ground contamination further. When the data themselves 
are analyzed as a whole, rather than separately for each 
modulation, the foreground analysis will also be done on 
the data set as a whole. This will take into account the 
cross-modulation correlations in the data. 

8. CONCLUSIONS 

The PyV experiment fully samples 598 deg^ of the mi- 
crowave sky and constrains the CMB angular power spec- 
trum in the angular scale range AQ <^ I <^ 260. The mea- 
surements pass internal consistency checks, show little con- 
tamination from foreground radiation, and are consistent 
with previous Python and COBE results. The observed 
angular power increases from larger to smaller angular 
scales, with a sharply increasing slope starting at I ~ 150. 
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Fig. 3. — Flat band power, {1(1 + l)C'i/2n)^'^ , vs. multipole I for all of the modulations. The detections have la error bars and the upper 
limit has 2cr error bars. The unapodized cosine modulation is plotted with an open square and the apodized cosine modulations are plotted 
with diamonds. The error bars include statistical uncertainties only and do not include uncertainties in the calibration or beam size. The I 
range of each modulation is determined by the half-maximum points of (W;)^". Low I values correspond to large angular scales and high I 
values correspond to small angular scales. CMB power is clearly rising from low to high / up to the sensitivity cutoff of PyV. 



